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METHOD FOR NEAR OPTIMAL JOINT CHANNEL EST IMATION 
AND DATA DETECTION FOR COFDM SYSTEMS 



Related Applications 

This application claims the benefit of priority of Provisional Application No. 
60/241,195, filed on October 19, 2000, and of Provisional Application No. 60/180,799, 
filed on February 7, 2000, and is related by subject matter to U.S. Patent application No. 
09/771,802 6.477.210, entitled "System for Near Optimal Joint Channel Estimation and 
Data Detection for COFDM Systems," by the inventors of the present application filed 
concurrently herewith. 

Field of Invention 

The invention relates generally to communications and particularly to a method 
and apparatus for near optimal joint channel estimation and data detection to improve 
channel tracking and, thus, improve link robustness. 



Background of the Invention 

The rapid growth in the use of the Internet and the increasing interest in portable 
computing devices have triggered the desire for high-speed wireless data services. One 
of the more promising candidates for achieving high data rate transmission in a mobile 
environment is Orthogonal Frequency Division Multiplexing (OFDM), which divides the 
wide signal bandwidth into many narrow-band subchannels, which are transmitted in 
parallel. Each subchannel is typically chosen narrow enough to eliminate the effects of 
delay spread. Coded OFDM (COFDM) systems, which combine both OFDM and 
channel coding techniques, are able to improve the performance further by taking 
advantage of frequency diversity of the channel. 

Though both differential and coherent demodulation can be applied in a COFDM 
system, the latter leads to a performance gain of 3 to 4 dB in signal-to-noise ratio (SNR) 
with accurate channel estimation. Channel estimation techniques realized by a 
frequency-domain filter using Fast Fourier Transform (FFT), followed by time-domain 
filters for a COFDM system with Reed-Solomon (RS) coding have been proposed. These 
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channel estimation techniques, while good, did not provide the near optimal channel 
estimation required for data-decoding with improved channel tracking capability for 
reliable link performance even under high user mobility and/or high RF carrier 
frequency. 

[ [Summary of the Invention] ] 

The physical layer configuration of the system for near optimal joint channel 
estimation and data detection for COFDM systems is depicted in FIG. 1 A. An exemplary 
transmitter 105 is shown on top and an exemplary receiver 140 is shown on the bottom of 
FIG. 1 A. At transmitter 105, a data stream is accepted by a convolutional encoder 1 10, 
which convolutionally encodes the data stream. The encoded data stream may then 
optionally be forwarded to an interleaver 1 1 5 for interleaving. If the encoded data stream 
is interleaved in the transmitter 105 then the receiver 140 must correspondingly de- 
interleave the encoded data stream. After interleaving, a modulator 120, for example a 
QPSK modulator 120, modulates the encoded (and optionally interleaved) data stream, 
which is then forwarded to an inverse Fast Fourier transformer 125, to subject the 
modulated encoded (and optionally interleaved) data stream to inverse Fast Fourier 
transformation. The transformed modulated (and optionally interleaved) encoded data 
stream (signal) is then transmitted, in the present invention, over the air via RF unit 130 
and antenna 135. 

Correspondingly, receiver 140 accepts multicarrier transmitted signals (data 
streams) via antennas 145 and RF units 150 and subjects the received multicarrier signals 
to Fast Fourier transformation using Fast Fourier transformers 155. These transformed 
signals are concurrently fed into channel estimator 165 and demodulators 160, for 
example QPSK demodulators. The demodulated transformed signals are combined in 
maximum ratio combiner 170. The combined demodulated transformed are then 
optionally de z interleaved using de-interleaver 175. The combined demodulated 
transformed (and optionally de^interleaved) signal is then decoded using Viterbi decoder 
180. The decoded combined demodulated transformed (and optionally de-interleaved) 
signal is then fed back into channel estimator 165, which forwards channel estimations, 
which are added to the transformed signals that are forwarded to demodulators 160. 
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FIG. IB shows the baseband processing, in particular, the iterative nature of the 
receiver portion of the system disclosed herein for near optimal joint channel estimation 
and data detection for COFDM systems. Channel estimator 165 accepts transformed 
signal 190. Channel estimations 198 are fed back into channel estimator 165. Channel 
estimations 194 are fed into decoder 185, which comprises maximum ratio combiner 170 
(shown in FIG. 1 A), optional de-interleaver 175 (shown in FIG. 1A) and Viterbi decoder 
180 (shown in FIG. 1A). Channel estimations 195 are fed into decoder 185 via 
demodulator 160 (shown in FIG. 1 A, but not shown in FIG. IB for clarity and to 
highlight the iterative nature of the system), which demodulates the transformed signal 
using channel characteristics. Decoder i9$ 185 also accepts transformed signal 190. 
Decoder 185 outputs signal 192, which is fed back into channel estimator 165. 

An OFDM signal is divided into a number of subchannels. By way of example, 
an OFDM signal bandwidth is divided into 120 6.25-kHz subchannels with QPSK 
modulation on each subchannel. At the receiver, the demodulated signals from two 
receiving branches are combined using maximal ratio combining and then decoded. With 
a symbol period of 200 ju s (including a 40- // s guard interval) and 1/2-rate coding, a 
maximum information rate of 600 kbps can be achieved in a 750-kHz bandwidth (about 
800 kHz including guard bands). The information rate is calculated by dividing the 120 
subchannels (tones) by the 200 n s period to obtain 600 kbps. 

For purposes of example for the present invention, 1/2-rate convolutional codes 
(CC) are considered. The results with 1/2-rate Reed-Solomon (RS) code based on 
Galois-Field (64) (GF(64)) are compared. The size of a code word is the same as that of 
an OFDM block (an OFDM symbol of 200 /u s and 120 subchannels). To achieve coding 
gain with inherent frequency diversity in OFDM, a simple interleaving scheme is applied. 
For both RS and CC cases, the first 120 bits of a code word are assigned to the in-phase 
component and the rest to the quadrature component. To gain additional randomness 
within a code word for the CC case, each 120-bit group is interleaved over subchannels 
by an 1 1-by-l 1 block interleaver (without the last bit). 

In the simulations, the wireless channel, as a Rayleigh-fading channel, with a two- 
ray multipath delay profile is modeled. Good performance for impulse separation as high 



5 



Chuang 2000-0071 



as 40- n s can be achieved; a 5- jus impulse separation in the numerical results is 
considered. 

[Oil] For the performance with respect to channel variations, maximum Doppler 

frequency up to 200 Hz, which is reasonable for most vehicular speeds, for a possible RF 
carrier frequency around 2 GHz is considered. To demonstrate the advantage of the 
proposed joint detection methods, results at a maximum Doppler frequency as high as 
500 Hz corresponding to a scenario in which the wireless system uses a higher carrier 
frequency, e.g. 5 GHz are presented. 

[012] In the medium access control (MAC) layer, a frequency reuse is considered with 

dynamic resource management, e.g., Dynamic Packet Assignment (DP A), to achieve 
high spectral efficiency for packet data access. 

[013] A simple analysis to highlight the ideal or optimal joint channel estimation and 

maximum likelihood (ML) decoding scheme indicated in FIG. 1 A for the case of M = 2 
receiving antennas is now presented. 

[014] At a diversity receiver, the signal from the m th antenna at the k th subchannel 

and the n th block can be expressed as 

where a t , h b and w m „ k are the transmit signal, channel response and additive 

ft AC * fTl yft yK ffl j" yft 

Gaussian noise, respectively. 
[015] For convolutional codes, because the size of a code word is the same as that of the 

OFDM block, (1) can be rewritten as 

V«= H *,A,+ W m ,«> (2) 
where, if there are K f subchannels, H W>M =rftog(^ nJ A^v-'A,!..^)' c « isthe 
transmitted code word at time epoch n , and the rest of the vectors are similarly defined. 
[016] Assume that the number of code words is N , we introduce the following 

notations, 

_r T T o r i r 

H m = diagiU,,, ^,...,^) , (3) 

-r r T x T ~\ T 
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[017] At the receiver, the objective is to solve a maximum likelihood (ML) problem 

c = argmin[minXl|x m -H m c|| 2 ], (4) 

with a constraint on channel response 

KHJ = 0, (5) 
where Z() is a constraint function. In a wireless environment, this constraint can be 
simplified to be 

EB„, / d(H m ,„_ / ) = 0, (6) 

where the length of the channel memory is K m OFDM symbol durations, B„ , are 
coefficients determined by the correlation between channel responses at the time epochs 
n and n-1, which is a function of the Doppler spectrum of the channel, and d() is a vector 

function defined by d(H OT „) = [h m „,vK^-^K^,K f f • 
[018] The optimal solution of this ML problem can be obtained by exhaustive search. It 

requires solving the mean square error (MSE) 

MSE(c) = min]>> m -H m c|| 2 , (7) 

for any possible c with the channel constraint (6). Then, 

c = argminMSE(c). ( 8 ) 

c 

[019] After obtaining MSE(c) , the corresponding channel estimate H m (c) can be 

found. Consequently, the optimal approach for estimating channel response requires the 
knowledge of the entire set of x and c . 

[020] Another observation from this ML receiver is that the channel estimation results 

H m is not a direct output of the detection process and hence, channel estimation which 
calculates H m explicitly may not be necessary in theory. However, for other required 
parameter estimation, such as timing and frequency synchronization, a known data 
sequence is usually transmitted in the beginning of a group of OFDM blocks. This 
known data sequence, also called a synch word or a unique word, can be used as a 
training sequence in (7) to obtain initial channel estimate explicitly without resorting to 
blind detection. This initial channel characteristic is helpful for solving this ML problem 
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with better numerical stability and tracking property. This initial channel estimation can 
be easily solved in the frequency domain by first taking FFT as shown in FIG. 1 A. 

[021] One related method and system is U.S. Patent N o. 6.327.314 Serial Numb e r 

09/089,862, METHOD AND APPARATUS FOR CHANNEL ESTIMATION FOR 
MULTICARRIER SYSTEMS, which was filed June 3, 1998, and is commonly held and 
incorporated herein by reference. The near optimal joint channel estimation and data 
detection method and system of the present invention was born from the research that 
resulted in that application and the subsequent determination that improvements could be 
made in the channel estimation. 

[022] The sub-optimal approach of the related system and method is now outlined. 

Because of the formidable complexity of the optimal ML receiver, some sub-optimal 
solutions are widely used in practice. The related sub-optimal solution is to divide the 
ML problem into two parts, channel estimation and coherent decoding. Then, the 
problem can be solved by iteratively estimating channel and decoding in the forward 
direction (in time). 

[023] At a time instant n , given a channel estimate H m „ , initially obtained by using 

training sequence (in the frequency domain dividing the transfer function of the received 
signal by the transfer function of the known data), the maximum likelihood (ML) 
problem 

c„ = argmin^||x m n -H m „c„ || 2 (9) 

can be solved. Then, the reference for channel estimation is 

H mjn =argmin2l|x m ,„-H m)n cJ| 2 - ( 10 > 

Finally, the estimate for the time instant n + 1 is obtained by solving a linear constrained 
equation (6). Considering a stationary channel with fixed maximum Doppler frequency, 
the coefficients B„ , are independent of n , and can be written as B, , which are used as 
the coefficients in the FIR filter to track channel variations. Consequently, a 
simplification of (6) with only previous references for prediction-type estimation leads to 

ZB / d(H m>n+l . / )-d(H m ,„ +1 ) = 0, (11) 

/=i 
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where M, is the number of taps of an FIR filter and B, are preset coefficients designed 
to achieve the minimum mean square error (MMSE) of estimation. This MMSE 
estimator can be realized by a frequency-domain filter using the Fast Fourier Transform 
(FFT), followed by an M, -tap time-domain filter, 

B/^F-'B^, (12) 
where b, is the time-domain filter coefficient, F is the FFT matrix, B f is a diagonal 
matrix, and F" 1 B / F is the frequency-domain filter. 

[024] The MMSE filter coefficients b, and B f were derived assuming c„ = c„ for a 

given set of Doppler frequency and delay spread. It is shown that this estimator is robust 
regardless of frequency or time mismatches. With a low Doppler frequency, it has been 
shown, that a 5-tap ( M, =5) estimator can successfully predict the channel. 

[025] To obtain accurate initial channel estimation, a training OFDM block is sent at the 

beginning of a transmission, in which c, is known to obtain H m , . The channel 
parameters for the new time epochs and the unknown code words can be successively 
solved in the forward direction (time advance). 

[026] The assumption c„ = c„ cannot be always guaranteed; an incorrectly detected 

code word introduces wrong channel estimation and hence, can cause a wrong detection 
of the successive code words. This error propagation is the dominant impairment of the 
link performance at high SNR. In order to alleviate the problem of error propagation, a 
training OFDM block is inserted every N, block. In the simulations presented herein, N, 

= W is considered. 

[027] In FIG. 2, the performance of the coherent reception with RS code and 

convolutional codes for 40-Hz maximum Doppler frequency is shown. The 
convolutional codes are shown with different constraint length (K) ranging from 3 to 9 as 
dashed lines . The performance of convolutional codes is substantially better than that of 
the RS code, which is shown as a solid line with cross lines. In order to achieve a Word 
Error Rate (WER) of 10' 2 , the K=9 CC needs 4 dB lower SNR than the RS code. In the 
use of WER in the present invention, it is assumed that a word is a codeword. Moreover, 
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the performance of the K=9 CC with channel estimation is very close to the one with the 
ideal channel information shown as a solid line . 

[028] In FIG. 3 , the performance at 200-Hz Doppler is shown. In comparison with the 

RS code, the CC's are still superior although the degradation with respect to the idealized 
case is higher due to poorer channel tracking. In fact, an error floor at the high SNR 
region exists due to tracking errors. Once again the RS coded signal is indicated as a 
solid line with cross lines. The convolutionallv coded signals are in dicated bv dashed 
lines and ideal channel information is indicated bv a solid line. 

[029] In FIG. 4, the performance of the K=9 CC with different maximum Doppler 

frequencies is shown. With a low Doppler frequency, a 5-tap (M, =5) estimator used 
here can successfully predict the channel and the performance is very close to that with 
idealized channel estimation. However, when the fading is relatively fast, it is difficult to 
estimate the channel correctly and the WER floors on the order of 10" 3 can be clearly 
observed for at a maximum Doppler frequency of 200 Hz. That is, it was found that the 
original method works well in slow fading but degrades significantly in fast fading. Once 
again the ideal Doppler frequency is indicated bv a solid line. A Doppler frequency of 
200 Hz is indicated bv a dashed line. A Doppler frequen cy of 175 Hz is indicated by a 
dashed line with small circles. A Doppler frequency of 15 0 Hz is indicated by a dashed 
line with triangles and a Doppler frequency of 125 Hz is indicate d bv a dashed line with 
cross lines . 

[030] FIG. IB depicts a block diagram for a baseband receiver. It is simplified to the 

extent that only channel estimator and decoder portions of the receiver are depicted. The 
decoder further comprises a plurality of demodulators, as well as a maximum ratio 
combiner, an optional deinterleaver (to match an optional interleaver in the transmitter) 
and a Viterbi decoder, all depicted in FIG. 1 A. However, FIG. IB depicts the structure 
and corresponding connections of the receiver depicted in FIG. 1 A. Note, all x, c , H, 
H are complex. Although it may seem strange to have a complex codeword, e , it is very 
natural to treat coding and modulation as a whole like spatial-temporal coding as 
discussed in "Space-time codes for high data rate wireless communication: performance 
criterion and code construction", by V. Tarokh, N. Seshadri, and A. Calderbank, 
published in IEEE Trans. Info. Theory, vol. 44, no. 2, pp. 744-765, Mar. 1998 or coded 
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modulation "Channel coding with multilevel/phase signals," by G. Ungerboeck, 
published in IEEE Trans. Info. Theory, vol. IT-28, no. 1, pp. 55-67, Jan. 1982. 

[031] At a time instant n, the channel estimator unit has two tasks. One is to produce 

the channel estimates of current time instant, H, „ , H 2n with its input xi,„, X2,„, ... and 
feedback H, „. M( ,H 2 n _ M( ,...,H /n . 1 ,H^„. 1 . by equation (1 1). The other task is to produce 
references H, n ,H 2 n for estimate processing in the next time instant by equation (10) 
when c„ is available. At a time instant n, when H ln ,H 2n ,... are available, the decoder 
unit produces decoded information c n by (9). 

[032] The flowchart of this method is shown in FIG. 1C. The related channel 

estimation method is first initialized at step 165-1. Transmitted signals are received at 
step 165-2 . A determination is then made at step 165-3 as to whether the received block 
is a training block. If the received block is a training block then c n is known and 

H m ,„ =argmin£||x m „ -H m n c„ || 2 
is calculated at step 165-5 . This is a reference for the channel estimation 

ZBAH m ,„ +w )-d(H m , n+1 ) = 0 

/=i 

which is calculated next at step 165-6 . The block number is incremented at step 165-7 
and a determination is made if the end of the frame has been reached at step 165-8. If the 
current block is not a training block then c„ is decoded at step 165-4 
c„ = argmin2^l|x m>n -H mi „c„ || 

C n m 

is calculated before calculating the reference and channel estimation. 
[033] It should be clear from the foregoing that there is room for improvement between 

the prior art and optimal joint channel estimation and data detection. An object, 
therefore, of the present invention is to improve the joint channel estimation. This will 
have the effect of reducing the impact of noise as well as reducing decoding errors. Thus, 
overall system performance will be improved. 
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[034] It is a further object of the present invention to provide a method and system that 

are robust even in light of a mismatch between Finite Impulse Response (FIR) 
coefficients and the true channel. 

[035] Another object of the present invention is to provide a method and system with 

improved channel tracking capability, resulting in reliable link performance even under 
high user mobility and/or high RF carrier frequency. With improved link performance, 
data rates significantly higher than currently available (or even than third generation 
systems in planning) can be offered to subscribers. 

[036] All of the above objects can be achieved nearly optimally even in rapid dispersive 

fading channels. 



Summary of the Invention 

An advance in the art is achieved bv realizing that an i terative approach for 
developing a channel estimation which employs current information, e.g. decoded d ata 
and newly received signal, yields a better channel estimation . At each time epoch n_, the 
improved receiver does not proceed to the time epoch n + ]_ immediately. Instead, the 
receiver uses the newly generated channel information for th e next time epoch, n+1, 
which is derived based on the current information, to substitute for the current channel 
estimate and repeats the calculation for the current channel est imation and the decoding 
of the transmitted symbols. In other words, the same procedu res are performed twice. 
Thus, an OFDM receiver of the present invention estimates the chann el bv receiving a 
multicarrier s i gnal, transforming the received signals bv applying Fast Fourier 
transformations, and estimates the channel characte ristics of the multicarrier channel 
using iterative processing, and in the process also decoding the transform multicarrier 
signal. 



Brief Description of the Drawings 
[038] FIG. 1 A is the physical configuration of a system for transmitting and receiving 

signals. 

[039] FIG. IB is a block diagram for a baseband receiver at time instant n. 

[040] FIG. 1C is a flowchart for a related sub-optimal channel estimation method. 
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[041] FIG. 2 is a comparison of RS and convolutional codes, where the maximum 

Doppler frequency is 40 Hz. 
[042] FIG. 3 is a comparison of RS and convolutional codes, where the maximum 

Doppler frequency is 200 Hz. 
[043] FIG. 4 depicts the performance of the K=9 CC at different maximum Doppler 

frequencies. 

[044] FIG. 5 depicts the performance of the iterative approach at 200 Hz maximum 

Doppler frequencies, where K=9. 
[045] FIG. 6A is a flowchart for the iterative processing. 

[046] FIG. 6B is a flowchart for the iterative backward processing. 

[047] FIG. 7 shows the performance of iterative backward-processing approach at 200 

Hz maximum Doppler frequency, where K=9. 
[048] FIG. 8 shows the performance of iterative backward-processing approach at 

different maximum Doppler frequencies, where K=9. 
[049] FIG. 9 is a comparison of systems having different time-domain FIR estimators, 

where K=9. 

[050] FIG. 10 shows the average retransmission probability as a function of occupancy 

per sector. Eight slots and three RF carriers (< 2.5 MHz) are reused in every base station, 
each with three sectors, using Dynamic Packet Assignment (DP A). In the original (non- 
iterative) method, K=5, 40 Hz Doppler or K=9, 125 Hz Doppler. 

[051] FIG. 1 1 shows the average delay of delivered packets as a function of throughput 

per base station, each with three sectors. Eight slots and three RF carriers (< 2.5 MHz) 
are reused in everywhere using DP A. In the original (non-iterative) method, K=5, 40 Hz 
Doppler or K=9, 125 Hz Doppler. 

[052] FIG. 1 2 depicts the average retransmission probability as a function of occupancy 

per sector. Eight slots and three RF carriers (< 2.5 MHz) are reused in every base station, 
each with three sectors, using DP A, where K=9, 200 Hz Doppler. 

[053] FIG. 13 shows the average delay of delivered packets as a function of throughput 

per base station, each with three sectors. Eight slots and three RF carriers (< 2.5 MHz) 
are reused in everywhere using DP A, where K=9, 200 Hz Doppler. 
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Detailed Description of the Preferred Em bodiments 
[054] A method and system using two channel estimation schemes, which improve 

channel-tracking capability based on simplification of the ideal or optimal detector are 
now presented. 

[055] The present invention uses a COFDM system with convolutional coding (CC). 

Moreover, the present invention uses two iterative processing techniques using current 
and future tentative data decisions in the backward direction for channel estimation and 
the final decoding. This permits joint channel estimation and data decoding with 
improved channel-tracking capability, resulting in reliable link performance even under 
high user mobility and/or high RF carrier frequency. System level performance such as 
retransmission probability and packet delay in a system with aggressive frequency reuse 
using dynamic packet assignment (DP A) is also presented. 

[056] Since the wireless channel is correlated in time, the optimal joint detection 

requires processing the received signals and the decoded data in the past, current and 
future epochs. From (10) and (1 1), it can be seen that the joint detector in the related is 
sub-optimal because it attempts to predict the channel response at a time instant by only 
using received signals and decoded information in the past. In order to improve channel 
estimation by applying the current information, e.g., decoded data and newly received 
signal, an iterative solution can be applied. 

[057] At each time epoch n , after performing (9), (1 0) and (1 1 ), the present invention 

does not proceed to the time epoch n + 1 immediately. Instead, the present invention 
uses the newly generated H m n+1 , which is derived based on the current information, to 
substitute for the current channel estimate H m n and the computations in (9), (10) and 
(1 1) are repeated. In other words, the same procedures are performed twice. This 
additional iteration allows channel estimation to be performed based on the current 
information. 

[058] In FIG. 5, the improvement of the system performance is shown. At each time 

instant n , if the perfect past channel were known, i.e., [H m ]) „H m „. i ] = [H m i„,H m „_,] in 
(11), the new iterative processing can perform within .3 dB from the case in which 
channel is known. On the other hand, even if the perfect past channel information were 
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available, the system performance would still be far from the optimal one with an 
irreducible error floor in the related sub-optimal solution. The SNR requirement of the 
new iterative approach at a WER of 1 0' 2 is 1 .2 dB lower than that of the original one. 
System performance for the oripinal signal is in d i cate d by a dashed line. Using an 
iterative ap proach is indicated bv a da shed line with triangles. Using the original 
ap proach and having per fect past chann e l information is indicated by a dashed line with 
circles T Ising an iterative ap proach a n d having perfect past channel information is 
indicated hv a dashed l ine with cross li n es. The ideal estimate is indicated by a solid line . 

[059] More than two iterations were attempted but the improvement was limited. 

Therefore, it was determined that another possible improvement probably should come 
from utilizing/wfure information as described next. 

[060] Following the insights of the optimal approach, by taking into account both the 

current and the future information another step is performed. To take advantage of the 
future information in the channel estimation process, the FIR filter shown in (1 1) can be 
processed in the time-reversed fashion, i.e., 

^B / d(H m , n+M )-d(H m , n .,) = 0. 03) 

/=i 

Hence, the iterative approach is performed backwards. Here, assume that the iterative 
approach mentioned in the last sub-section has been already carried out from time instant 
1 to N . Therefore, x m and c are available with which to perform the iterative approach 
backward at a later time for tentatively detected previous signals; this method is thus 
referred to as iterative backward-processing approach. In order to process backward, 
H is stored when the decoding result is correct (assuming an error detection 
mechanism is available) and c and x m when a decoding error occurs at a time instant n . 
[061] In general, there are three choices to perform the backward processing. The first 

choice is to process backward after all N OFDM-block data are processed. In reality, it 
is not a good choice because the channel memory length K m is not infinite and the delay 
of this approach may not be acceptable if N is large. A feasible approach is to process 
backward after having M, (FIR tap length) consecutive correct code words (OFDM 
blocks). This ensures that the backward processing will start from a reliable channel 
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estimate. The last choice is to process backward starting at any other points, which is 
determined by memory and delay requirements. Performance degradation is expected. 
However, if the backward processing is started from a known training OFDM block, 
lower degradation can be achieved. With the last two choices, it is only necessary to 
store M N (« N in general) OFDM-block data in the memory. 
[062] The iterative processing and iterative backward processing approaches share the 

similar baseband processing block diagram (FIG. 1 A) with the related/original approach. 
The only difference is that the feedback is H, „. M+1 ,H 2 n . M+ ,,...,H, „ ,H 2 „ and 
G u H. , H n , , respectively, when performing iterative processing and 

iterative backward processing. Their corresponding flowcharts are shown in FIG. 6A and 
FIG. 6B. 

[063] The iterative channel estimation method depicted in FIG. 6A is first initialized at 

step 605 . The iterative processing for estimating channel ch aracteristics is performed by 
using the system depicted in Figs. 1 A and IB and a s described above. Transmitted 
signals are received at step 610 . A determination is then made as to whether the received 
block is a training block at step 615 . If the received block is a training block then c„ is 

known and 

H m „ = argmin£||x m> „ -H m „c„ || 2 

is railrailate d. which is a tentative reference signal, bv first tentatively decod ing the block 
of the received multicarrie r si gnal at step 620 . This is a reference for tho channel 
estimation. The tentative reference signal is then use d to generate a tentative estimation 
for the channel at step 625 given bv the equation 

E B / d (H m ,« + .-/)- d ( H m , + .) = 0 



/=i 



which is calculated next. The tentative reference signal is then used to generate a 
tentative estimation for the channel at step 625. the block number is incremented alstep. 
630, and a determination is made if the end of the frame has been reached a t step 635 . If 
the end of the frame has not been reached then another b lock of the received multicarrier 
i gnal is accepted for processing at step 610 . If the current block is not a training block 
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then 

c„ = argmin]£||x ffl „ -H m „c„ || 2 

is calculate d, which is a reference signal, by first decoding the block of the received 
multicarrier s i gnal at step 640. the matrix Then the reference 
H m n = argmin£||x mn -H m>n c„ || 2 

H m,n m 

is calculated followed by calculating at step 645, and it is then used to generate an 
estimation of channel characteristics given by 

ZB ; d(H m ,„ +1 _ / )-d(H m ,„ +1 ) = 0 

/=i 

thn nhnnnnl nnrimitinn at ste p 650. The block of the received multicarrier signa l is then 
re-decoded using the estimation of the channel char acteristics at step 655. The method 
then proceeds to the atop after determining that the current block ia a training block step 
620. This effectively repeats the reference and ch annel estimation. The reference and 
channel estimation are repeated in order to improve the calculations with the tentative 
reference and channel estimation calculations. 

The iterative backward channel estimation method depicted in FIG. 6B is first 
initialized at step 660 . The iterative backward p rocessing for estimating channel 
characteristics is performed bv using the system de picted bv Figs. 1 A and IB and as 
described above. Transmitted signals are received at step 665 . A determination is then 
made as to whether the received block is correct at step 670 . If the received block is 
correct then c„ is known and 

H m n = argmin£||x ffl „ -H m n c„ || 2 

is calculate d, which is a tentative reference signal, bv firs t tentatively decoding the block 
of the received multicarrier signal at step 675. This ten tative reference signal is then used 
to generate a tentative estimation for the channel at st e p 680 given by This is a reference 
for the channel estimation 

/=i 

which io calculated next . The block number is decremented at step 685 and a 
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determination is made if the beginning of the frame has been reached at step 690 . Ifthe 
beginning of the frame has not been rea c hed then another block of the received 
mnlticarrier signal is accepted fo r processing at step 665. If the current block is not 
correct block then £„ 

c„ = argmin^llx^ -H.^c, || 2 

is ralrniate ri, which is a reference signal, bv fir st decoding the block of the received 
multicarrier signal at step 692 . Then-the This reference giyenby 
H m „ = argmin^||x m> „ -H m „c„ || 2 

is calculated followed by calculating at step 694 is then used to generate an estimation of 
channel characteristics at step 696 given by 

lB / d(H m , n+/ . l )-d(H m , n - 1 ) = 0. 

/=i 

tho channel estimation. The current block io then rcdocodod. The method thon proce e ds 
to the otop after determining that the current block is correct. The block of the received 
multicarrier sipnal is then re-decoded using t he estimation of the channel characteristics 
at ste p 698. and control is returned to step 675. This effectively repeats the reference and 
channel estimation. The reference and channel estimation are repeated in order to 
improve the calculations with the tentative reference and channel estimation calculations. 

In FIG. 7, the performance of this approach (iterative backward processing) is 
shown for 200 Hz maximum Doppler frequency. In the simulation, the maximum M N is 
set to be 200 and this corresponds to 40 ms. However, due to the low error probability in 
the high SNR region, much shorter storage is required. For instance, the maximum M N 
required is about 50 at the 5 dB SNR. It is found that nearly optimal performance is 
achieved with iterative backward processing. System performance using the original 
ap proach is indicated bv a dashed line. Sys t em performance using the iterative approach 
is indicated bv a dashed line with triangl es. System performance using the iterative 
backward processing approach of FIG. 6B i s indicated bv a dashed line with cross lines. 
The ideal estimate is indicated bv a solid line. 
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[066] In FIG. 8, the performance of this approach with different maximum Doppler 

frequencies is shown. The system still performs well in an environment with maximum 
Doppler frequency as high as 400 Hz. Once again the ideal estimate is indicated by a 
solid line. Syste m performance at 500 Hz usin p iterative backward processing is 
indicated hv a dashed line with crosses (or VsV System performance at 400 Hz using 
iterative backward proce ssing is indicated hv a dashed line with small circles. System 
performance at 300 Hz using iterative b ackward processing is indicated by a dashed line 
with trian ples. System performance at 2 00 Hz using iterative backward processing is 
ir.Hir.ated bv a dashed line with cross lin es. System performance of the original approach 
at 200 Hz is indicated bv a dashed line. 
[067] By applying the simple iterative estimation, the impact of noise is reduced by an 

additional round of filtering with newly available information. However, it may not 
eliminate the impact of decoding errors. Using iterative backward-processing, there is a 
chance to re-estimate the channel with fewer decoding errors and hence, better 
performance is achieved. 
[068] It should be noted that the Finite impulse Response (FIR) estimator coefficients 

B, in the related sub-optimal method were optimized for the scenario without any 
iterative process. It should be noted, however, the detector architecture introduced here 
can be combined with any filter design for better tracking. 
[069] To show the robustness of the iterative approaches, a simple averaging 5-tap FIR 

filter (b, = . 2) for the time domain filtering is now considered. As shown in FIG. 9, 
even with these simple FIR coefficients, the method of the present invention still 
outperforms the original method with the FIR that was optimized for a particular set of 
maximum Doppler frequency and delay spread. Therefore, the iterative backward- 
processing approach of the present invention is relatively robust against the mismatch 
between the FIR coefficients and the true channel. The ideal estimate is indicated by a 
solid line. System performance usi n g the original a ppr oach and an optimal 5-tap FIR at 
700 H7. is indicate bv a dashed and dotted line. System performance using the original 
ap proach and an optimal: 5-tan FIR at 40 Hz is indicated bv a dashed line. System 
performance using iterative backwa rd processing and an optimal 5-tap FIR at 200 Hz is 
indicted bv a dashed line with crosses ( or "x"sV System performance using iterative 
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W.kward processing and an optimal 5-tap FIR at 40 Hz is indicated by a s olid line with 
crosses for"x"'s). 

[070] Consider a MAC layer configuration to characterize the system-level performance 

under frequency reuse using the improved detection methods. A simulation system of 36 
base stations arranged in a hexagonal pattern is used, each having 3 sectors, with a 20-dB 
front-to-back ratio and idealized antenna pattern. The same channel can be used 
everywhere, even in different sectors of the same base station, as long as the Symbol to 
Interference Ratio (SIR) in the DPA admission process exceeds 7 dB. The simulation of 
the simultaneous use of the same spectrum by users in different cells results in 
interference between cells. Interference is treated as if it behaves like noise. 

[071] A channel is defined to be a combination of time slot and RF carrier, each 

consisting of 120 subchannels described previously. By way of example, the results for 
the case with 3 RF carriers and 8 time slots are now presented. This occupies a total 
spectrum of less than 2.5 MHz, including guard bands and other overhead. Each time 
slot consists of 10 OFDM blocks, i.e., 2 msec. One of these OFDM blocks is used for 
training, as discussed previously, while an additional block can be allocated for guard 
time between time slots. In addition, assume that a control slot of duration 4 msec is 
inserted in the beginning of every frame of 8 traffic slots to enable paging, assignment 
and pilot transmission that are required for the DPA process. With this conservative 
assumption of overhead in time and frequency domains, 48 kb/s (8 OFDM blocks or 960 
data bits transmitted in 20 msec) can be delivered using each time slots. Once paged, a 
mobile station (MS) measures the pilot signals to determine the desired traffic slots and 
reports the list back to the base station (BS). The BS then assigns traffic channel(s) and 
informs the MS this assignment for traffic packet delivery. 
[072] Based on the downlink frame structure, four adjacent BS's form a reuse group and 

they take turns performing the DPA procedure once every 4 frames. The time-reuse 
groups in the entire service area are pre-planned in a fixed and repeated pattern. 
[073] For the propagation model, the average received power decreases with distance d 

as d" 4 and the large-scale shadow-fading distribution is log-normal with a standard 
deviation of 10 dB. A data-service traffic model, based on wide-area network traffic 
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statistics, which exhibit a "self-similar" property when aggregating multiple sources, was 
used to generate packets. 

Automatic Retransmission reQuest (ARQ) is employed for retransmission when a 
packet is received in error. A packet in this case is theoretically 8 code words in each 
time slot, but the error probability is represented by using the WER curves. Since the 
error probability of the 8 code words in a time slot is highly correlated and additional 
coding is usually included for the entire packet, this approximation provides reasonable 
performance estimation for the MAC layer. If a packet cannot be successfully delivered 
in 3 seconds, which may be a result of traffic overload or excessive interference, it is 
dropped from the queue. The control messages are assumed to be error-free in the 
designated control slots. 

First consider K=5 and 40 Hz or K=9 and 125 Hz based on the related/original 
(non-iterative) method, both give similar link performance (see Figs. 2 and 5), for a 
comparison between RS and CC coding methods. FIG. 10 shows the average probability 
of packet retransmission, as a function of occupancy for all available (24) channels in 
each sector. This is a measure of QoS (quality of service) experienced by individual 
users. With a 3-6 % target retransmission probability, 15-50 % occupancy per radio in 
each sector is possible with this DPA scheme, depending on the use of coding schemes. 
Clearly, joint channel estimation and maximum likelihood detection of CC indicated by a 
dashed line introduced previously provides significant improvement over the case of RS 
codes, indicated by a solid line with cross lines which is also similar to the case of 
differential demodulation of the RS codes with 4 transmit antennas. Both results are 
significantly superior to the efficiency provided by current cellular systems, which are 
typically designed for voice communications with very conservative frequency reuse, 
about 4-7% spectrum occupancy in each sector. Data applications, permitting some 
retransmission delay, and improved link design, introduced here, allow much more 
aggressive and efficiency frequency reuse. 

FIG. 1 1 shows that 1-1 .5 Mb/s can be successfully delivered by each base station 
with an average delay on the order of 40-120 msec. This is a measure for system 
capacity. It indicates that OFDM link and DPA MAC combined enable a spectrally 
efficient (40%-60% b/s/Hz with a conservative assumption of overhead requirements) air 
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interface for broadband services, even for the macrocellular environment considered here. 
Adaptive modulation has not been considered in this study, and its use is expected to 
improve efficiency beyond 1 b/s/Hz per base station even under aggressive frequency 
reuse. The OFDM technology discussed herein can provide robust performance with 
peak-rates scalable with the available bandwidth. RS codes are indicated by a solid line 
with cross lines. CC codes are indicated bv a dashed line. 

Next, consider the case of high maximum Doppler frequency (200 kHz) and K=9 
(WER curves in FIG. 6) for comparison between sub-optimal ("original") detection 
method discussed earlier and near-optimal iterative backward-processing method 
presented. FIG. 12 shows that retransmission probability using the improved method can 
work well even under high maximum Doppler frequency. As a result, QoS can be 
improved even for high mobility users or when higher carrier frequency is employed. 
System performance using the original detection method i s indicated by a solid line with 
cross lines. System performance using the iterative backward processing approach is 
indicated by a dashed line. On the other hand, the capacity difference is relatively 
smaller, as shown by the delay-throughput curves in FIG. 13. System performance using 
the original a pproach is indicated bv a solid line with cross line s and system performance 
using the iterative backward processing is indicated b v a dashed line This is because both 
methods give very good radio link performance and the delay is dominated by sharing 
limited number of traffic slots, which is independent of the WER performance. In this 
case, better traffic resource management, such as improved admission control, could 
achieve capacity improvement. If higher Doppler frequency, e.g., 400 Hz, were 
encountered, the improved link performance introduced by the iterative backward- 
processing method, as shown in FIG. 8, would also result in system capacity 
enhancement. 

Returning to FIG. 1 A to show the detail of the decoder (DEC) unit and its relation 
with the Channel Estimator (CE) unit, with channel estimates H m „ , c„ can be calculated 
by maximum likelihood (ML) decoding (9), 

c„ = argmin^TlIx^ -H m „c„ || 2 . 
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[079] In FIG. 1 A, this decoding process is divided into 4 parts, e.g. QPSK 

Demodulator, a maximum ratio combiner (MRC), Deinterleaver and Viterbi decoder. 
How to separate these 4 units is now explained. 

[080] We have 

Zlk,-H m ,,cJ| 2 

m 

k m 

= E©lw II 2 - 2Re {<«A,^,*HI W,..* II 2 )- 



k m 



Note that J^IK, „ k || 2 makes no contribution to our minimization and due to binary 

k m 

convolutional code with QPSK modulation, || c n k \\ 2 is a constant. 

X^X'I Knk\\ 2 \\ c »k II 2 ) niakes no contribution to our minimization either. So we 



k m 

concentrate on 



k m 

r 



-z 



Re{c„,}- 



Zim{*L.«A, n ,*} 



Im{c„, t } 



The demodulator and the MRC calculate ^Ref^A,*,*} m<i Z Im {*».»A.»,*}- 

And the remainder is done by the deinterleaver and the decoder. 
[081] Conceptually, the MRC does only energy combining. So, we perform the 

following calculations: 

^Im^A,*,*} = X 1 " 1 ^.",* "BCw)} II II • 
Then, Re{x; n ,arg(^ nJt )} and Im{x; „ A arg(^ n t )} are actually QPSK demodulation. 

And the remainder is the MRC. 
[082] The radio-link performance of our COFDM system with the novel and non- 

obvious combination of two channel estimation schemes under fast fading validate the 
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design methodology and system for near optimal joint channel estimation and data 
detection. 

For a COFDM system, coherent detection can significantly improve radio link 
performance, but it requires a channel estimator to achieve the potential gain. The 
present invention applies a forward channel estimation and maximum likelihood 
decoding method for a COFDM system with convolutional codes and a simple 
interleaving scheme. It is shown that the performance of such a system significantly 
outperforms the one with Reed-Solomon codes. The performance under different fading 
rates was studied to understand the performance limits and areas for improvement. 
Simulations showed that this detection method is able to achieve near optimal 
performance in a wireless environment with maximum Doppler frequency as high as 100 
Hz. However, as fading rate increases further, irreducible error floor is introduced by the 
limitation in channel-tracking capability. A near-optimal detection method has been 
presented to improve the channel tracking performance, which is based on iterative 
processing and iterative backward processing of channel estimation and data decoding. 
Within practical values of interest, error floor can be eliminated even for very high fading 
rates, which is important for cases with high mobility or when higher carrier frequencies 
are considered in the future. Radio system performance based on the improved link 
techniques and a medium access control protocol using dynamic packet assignment was 
also studied. With radio link improvement, system throughput and delay can be 
significantly enhanced for better quality of service and spectrum efficiency. The system 
is also greatly simplified by allowing reuse factor of one without frequency planning. 

It should be clear from the foregoing that the objectives of the invention have 
been met. While particular embodiments of the present invention have been described 
and illustrated, it should be noted that the invention is not limited thereto since 
modifications may be made by persons skilled in the art. The present application 
contemplates any and all modifications that fall within the spirit and scope of the 
underlying invention disclosed and claimed herein. 



24 



